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ABSTRACT 

The emission line luminosity function of active galactic nuclei (AGN) is measured 
CN \ from about 3000 AGN included in the main galaxy sample of the Sloan Digital Sky 

I Survey within a redshift range of < z < 0.15. The Ha and [OIIIJA5007 luminosity 

functions for Seyferts cover luminosity range of 10^~^Lo in Hq and the shapes are well 
fit by broken power laws, without a turnover at fainter nuclear luminosities. Assuming 
C3 ■ a universal conversion from emission line strength to continuum luminosity, the inferred 

B band magnitude luminosity function is comparable both to the AGN luminosity func- 
tion of previous studies and to the low redshift quasar luminosity function derived from 
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O \ the 2dF redshift survey. The inferred AGN number density is approximately 1/5 of 

\ all galaxies and about 6 x 10^'^ of the total light of galaxies in the r-band comes from 

' the nuclear activity. The numbers of Seyfert Is and Seyfert 2s are comparable at low 

luminosity, while at high luminosity, Seyfert Is outnumber Seyfert 2s by a factor of 2-4. 
^ ' In making the luminosity function measurements, we assumed that the nuclear lumi- 

. nosity is independent of the host galaxy luminosity, an assumption we test a posteriori, 

and show to be consistent with the data. Given the relationship between black hole 
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mass and host galaxy bulge luminosity, the lack of correlation between nuclear and host 
luminosity suggests that the main variable that determines the AGN luminosity is the 
Eddington ratio, not the black hole mass. This appears to be different from luminous 
quasars, which are most likely to be shining near the Eddington limit. 

Subject headings: galaxies: active — galaxies: Seyfert — galaxies: starburst — galaxies: 
quasars: emission lines — surveys 

1. Introduction 

Active Galactic Nuclei (AGN), including high-luminosity quasars and low-luminosity Seyferts 
or LINERs, are among the most intriguing objects in the universe. The optical luminosity functions 
of AGN overall, as well as different types of AGN, hold important clues about the demographics of 
the AGN population, which in turn provide strong constraints on physical models and evolutionary 
theories of AGN. Many attempts have been made to derive the luminosity function of quasars 
(Schmidt & Green 1983; Marshah et al. 1983; Boyle et al. 1988, 1990, 2000; Hewett, Foltz & 
Chaffee 1993; Schmidt, Schneider & Gunn 1995; Hawkins & Veron 1995; Koehler et al. 1997; 
Goldschmidt & Miller 1998; La Franca & Cristiani 1997, 1998; Fan et al. 2001). The largest 
study to date is that of Groom et al. (2004), who evaluate the QSO luminosity function and its 
cosmological evolution from the combined 2dF QSO Redshift Survey and 6dF QSO Redshift Survey. 
Their luminosity function covers a luminosity range of Mb < — 22.5(with Hq = lOOkms"^ Mpc~^) 
and redshift range of 0.4 < z < 2.1. However, selecting AGN via their colors requires that the 
nuclear luminosity be at least comparable to the host galaxy luminosity to be detected. Therefore 
this method is biased against low luminosity AGN. 

To extend the AGN optical luminosity function to low luminosities, where the contribution of 
the host galaxy can equal or exceed that of the AGN, we must select AGN via their spectroscopic 
features in a galaxy redshift survey. Huchra & Burg (1992) selected 25 Seyfert Is and 23 Seyfert 2s 
from the CfA redshift survey (Huchra et al. 1983), and used these AGN to measure their luminosity 
functions. More recently, Ho et al. (1997a) carried out a careful and uniform spectroscopic survey 
of 486 nearby galaxies selected from the RSA catalog (Sandage & Tammann 1981) using very small 
apertures centered on the nuclei. From these galaxies they found 211 AGN, including 94 LINERs, 
65 transition objects and 52 Seyferts. Ulvestad & Ho (2001) obtained an optical AGN luminosity 
function from their sample. Vila-Vilaro (2000) used this sample as well as the starforming galaxy 
sample to estimate the shape of the luminosity functions for Seyferts, LINERs, transition objects 
and starforming galaxies. The luminosity functions obtained from these analyses extended to much 
lower AGN luminosities, Mb ~ —17, than that obtained from color-selected AGN samples, but they 
suffer from two disadvantages: 1) the AGN samples arc small, therefore the luminosity functions 
obtained from these samples are subject to large uncertainties, and 2) previous studies all used 
the B band magnitude of the entire galaxy instead of the nuclear luminosity as the variable in 
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their luminosity functions. As a result, the luminosity functions are severely contaminated by the 

AGN host galaxies. In order to overcome these disadvantages, large redshift surveys and careful 
evaluation of the luminosity functions arc needed. The Sloan Digital Sky Survey (SDSS), with its 
large number of high quality spectra available, provides us a great opportunity to finally fulfill this. 

The SDSS (York et al. 2000; Stoughton et al. 2002; Abazajian et al. 2003, 2004) is an imaging 
(Gunn et al. 1998) and spectroscopic survey that will eventually cover approximately one-quarter 
of the Celestial Sphere and collect spectra of ~ 10^ galaxies and 10^ quasars. Software pipelines 
automatically perform the data processing: astrometry (Pier et al. 2003); source identification, 
deblending and photometry (Lupton et al. 2001); photometric calibration (Fukugita et al. 1996; 
Smith et al. 2002); spectroscopic target selection (Eisenstein et al. 2001; Strauss et al. 2002; 
Richards et al. 2002); spectroscopic fiber placement (Blanton et al. 2003a) and spectroscopic data 
reduction. The SDSS spectra are of high quality with spectral resolution of about 2000 and typical 
signal to noise ratio of 16 per pixel for the galaxy sample. This allows construction of a uniform 
large sample of AGN identified from their spectroscopic features. 

Hao et al. (2005) (Paper I) described the selection of AGN from the SDSS spectroscopic 
data. We start from a low-redshift {z < 0.33) galaxy sample that is complete in r band Petrosian 
magnitude at 17.77 (Strauss et al. 2002) over 1151 square degrees (see Zakamska et al. 2003 for 

a selection of high redshift type II quasars from the SDSS). To properly measure the emission-line 
properties of these galaxies, we applied a stellar-subtraction procedure to each galaxy in the sample 
with a set of absorption templates developed from a Principal Component Analysis (PCA) of SDSS 
absorption-line galaxies. The broad-line AGN are identified by their broad Ha emission line width: 
they all satisfy Fuh Width at Half Maxima of Ha (FWHM( Ha))>l,200 km s'^, a natural break 
point in the distribution of Ha line widths. The narrow-line AGN are selected via their locations 
in the emission-line ratio diagnostic diagrams (Veilleux & Osterbrock 1987). Kewley et al. (2001) 
developed theoretical upper-limits for starburst galaxies in the diagrams and proposed them to be 
used as the criteria to separate narrow-line AGN from the usual star-forming galaxies. Kauffmann 
et al. (2003) plot the SDSS galaxies in the [0III]/H/3 ~ [Nil]/ Ha diagram and noticed that the 
galaxies naturally separate into two branches. They suggest classifying the galaxies in the second 
branch with higher [Nil] /Ha value as AGN. The separation point of the two branches lies well 
below Kcwlcy's separation line in the diagram. Therefore, the Kauffmann criterion selects many 
more narrow-line AGN than Kcwlcy's criteria. In Paper I, we applied both criteria and compiled 
an AGN sample containing 1317 broad-line AGN, 3074 narrow-line AGN via Kewley's separation 
and 10,700 narrow-line AGN via Kauffmann's criteria. Because they are basically low-luminosity 
AGN, in this paper we will sometimes call these AGN Seyferts and refer to the broad-line AGN 
sample as the Seyfert 1 sample, the narrow-line AGN sample selected via Kewley's separation as 
the Seyfert 2 sample and that selected via Kauffmann's criteria as the Seyfert 2* sample. 

To measure the luminosity function, we first need to evaluate the selection function; this is 
discussed in §2. We define our emission line AGN luminosity functions and describe our method 
to measure them in §3. In §4 we show the luminosity function results. In §5, we apply several a 
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posteriori checks on our luminosity function results and the assumptions used in deriving them. We 

further compare our results with the AGN luminosity functions obtained from previous studies in 
§6. In §7 we discuss issues related to the reddening correction and we summarize in §8. Throughout 
this paper, we will be using a cosmology with (^1^,0,^) = (0.3,0.7) and Hq = lOOkms"^ Mpc~^. 

2. Selection Function 

SDSS spectra are taken with a fixed 3" aperture, thus the observed spectra represent the sum 
of both the nuclear light and part or all of the host galaxy light. The continuum component of the 
low-luminosity AGN are difficult to separate from the host galaxy spectrum. Therefore we choose 
to use the emission-line luminosities to represent the nuclear luminosities and as variables in the 
AGN luminosity function. Emission lines, such as [GUI], Ha and [Oil] are good representatives of 
nuclear activity, except that in some galaxies the emission lines also include a contribution from 
star-forming activity. As explained in Paper I and further in Section 4, this contamination is not 
severe for AGN identified via Kewley's criteria, since they are galaxies that are dominated by AGN 
components. AGN selected via Kauffmann's criteria, however, might suffer this contamination since 
many of them are AGN + starburst galaxies. As demonstrated in Kauffmann ct al. (2003), the 
[QUI] emission lines for these galaxies, however, include little contamination from star formation. 
Thus, the [QUI] luminosity function derived from the narrow-line AGN sample using the Kauffmann 
et al. criteria will also be a good indicator of AGN activity. In this paper, we will obtain the Ha, 
[QUI] and [Oil] luminosity functions for samples obtained by the criteria both of Kewley et al. 
(2001) and Kauffmann et al. (2003). 

As described in Paper I and Section 1, the selection criteria for our AGN include both an overall 
galaxy magnitude cut and selection via certain spectral features. Thus the detection probability will 
depend not only on the emission line luminosity, but also the spectral features of the host galaxy, 
associated noise, and the overall galaxy luminosity, which are different for each AGN. Therefore we 
will use Monte-Carlo simulations that include all these factors to obtain the detection probability 
function for each AGN. 

The basic idea of the simulation is to evaluate the probability of classifying a given object in 
our sample as an AGN, as a function of the nuclear luminosity. We make a simplifying assumption 
to allow us to continue: we assume that the nuclear luminosity of an AGN (as represented by 
emission-line luminosity) is independent of its host galaxy luminosity. Therefore, as the nuclear 
luminosity changes in the simulation, the properties of the host galaxy and its associated noise are 
kept unchanged. This assumption is not obviously true and has significant physical implications of 
its own. We will show in §5 that our data are consistent with this assumption. 

Based on the selection criteria described above, there are several possible reasons that an AGN 
will not be identified as the nuclear luminosity is decreased. First, if the nuclear luminosity is too 
weak, the emission lines will be undetectable. Second, as the emission line strength is decreased, 
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the noise in the measurement of the emission hne ratios and the Ha emission Une width will 

increase. For AGN near the limit of the selection criteria, this noise will sometimes cause them to 
be no longer identified as AGN. Third, as the nuclear luminosity is decreased, the overall flux of the 
galaxy, which is the sum of contributions from the nucleus and the host galaxy, is decreased as well. 
When the total brightness reaches the galaxy target magnitude limit ( r band Petrosian magnitude 
equal to 17.77) the AGN will drop out of our sample. Therefore, the detection probability increases 
as the nuclear luminosity increases. The detailed behavior depends on the spectral features of 
individual spectra, which are different from object to object. Therefore our simulations are carried 
out on an object to object basis. 

If we denote the original spectrum of an AGN observed in a 3" aperture as /(A) and the 
emission line spectrum obtained after the continuum subtraction (Paper I) as /e(A), then we can 
build a pure emission- line spectrum /e(A) from /e(A) that contains only the emission lines relevant 
to the AGN identification by multiplying /e(A) by a window function: 



Here, the A^ are the central wavelengths corresponding to Ha, H/3, [OIII]A4959 [OIII]A5008, 
[Nn]AA6548,84, [Sn]AA6719, 32 and [OI]A6300. Ak = Sa^, where Uk is the Gaussian width of 
the line k. If two components are preferred in the fitting for Ha or H/?, we use the width of the 
broader component (Paper I). 

We use /e(A) to represent the nuclear spectrum and /a(A) = /(A) — /e(A) represents the host 
galaxy spectrum. During the simulation, we evaluate the corresponding detection probability while 
changing the nuclear luminosity, i.e., the flux of the emission line spectrum /e(A), by scaling by a 
constant c < 1. Based on our assumption of no correlation between AGN and galaxy luminosity, 
the host galaxy spectrum /a(A) and the associated noise should be kept unchanged. However, when 
/e(A) is decreased, the noise within the relevant emission line range is scaled down accordingly, even 
though the host galaxy, which may dominate the spectrum, is unchanged. Therefore to keep the 
noise within the emission line range unchanged, we generate a noise spectrum. The noise associated 
with each wavelength pixel is randomly generated following a Gaussian distribution. The amplitude 
of the noise is chosen to be comparable to the noise of the original spectrum at each pixel. The noise 
spectrum is multiplied by the emission-line window function of Equation 1, thus the host galaxy 
spectrum outside of the emission-line range is unchanged. If we denote the final noise spectrum 
as /noise (A), we can build a new spectrum with decreased nuclear luminosity and emission line 
strength as: 



The whole process of AGN identification, including stellar subtraction, is applied to this new 
spectrum to check if it still identified as an AGN, i.e., for broad-line AGN, if its Ha emission line 
has FWHM larger than 1,200 km s~^ and for narrow- line AGN, if its emission- line ratios satisfy 
the AGN emission line ratio criteria (Kewley et al. 2001; Kauffmann et al. 2003; Paper I). 




1 — \k\ < Ak for any k 

otherwise 



(1) 



fnew{^) — /a(A) + fnoiseW + C ■ /e(A) 
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For each specific nuclear luminosity tested, i.e., for each value of c, the whole process is re- 
peated ten times, each with a different randomly generated noise spectrum /noise (A). The detection 
probability associated with this nuclear luminosity and therefore this coefficient c is given by the 
number of times the new spectrum is identified as an AGN, divided by the total number of tests, 
namely ten. 

The detection probability as a function of nuclear luminosity therefore can be obtained by 
carrying out this simulation for a series of nuclear luminosities, i.e., a series of coefficients c. The 
values of c are chosen to allow luminosities between ~ lO^L© and ~ 1O^L0 in logarithmic steps of 

]^q0.034 

The simulation is done for each AGN in our sample, resulting in a detection probability function 
as a function of c. For most AGN, the detection probability function is a monotonic function of 
AGN luminosity, from zero at low nuclear luminosity, to unity when the nuclear luminosity is equal 
to or less than the original nuclear luminosity, as shown in Figure 1(a). For some other AGN 
like that in Figure 1(b), the simulation noise overwhelms the signal and the detection probability 
function is not monotonic, but fluctuates before it reaches unity. This often happens to AGN near 
the limits of the selection criteria, i.e., for narrow- line AGN, these are galaxies located close to 
the AGN/starforming galaxy separation lines in the emission-line ratio diagram; or for broad-line 
AGN, they have Ho; FWHM close to 1,200 km s~^. As the nuclear emission lines become weaker, 
the noise becomes more prominent and begins to affect the measured emission line properties. 

Some AGN have detection probability less than unity even when the simulated nuclear lumi- 
nosity has reached the original nuclear luminosity, as shown in Figure 1(c). One would imagine 
that there is a luminosity higher than that of the original AGN at which the detection probability 
would flnally reach unity. However, it is hard to correctly estimate the detection probabilities for 
nuclear luminosities larger than the original nuclear luminosity, because of the difficulty of cor- 
rectly estimating the true spectrum at higher S/N. As a result, we will simply set the detection 
probabilities at nuclear luminosities larger than the original one to be the value at the observed 
luminosity. 

The detection probability as a function of c discussed above is calculated purely based on 
detectability of spectral features. In addition, nuclear continuum in broad-line AGN contributes to 
the broad-band flux. Some AGN entered our sample because the nuclear continuum contribution 
pushed their r band magnitudes above the magnitude limit of 17.77 of our sample. To quantify this 
effect, we need to have some knowledge of the non-stellar nuclear continuum contribution which we 
calculate from the emission line intensities. 

If we assume that the emission lines we observed in our AGN spectra are all coming from the 
nucleus, then they are excited by the nuclear continuum light generated by the accretion. Therefore, 
the intensity of the emission lines and their continuum must be correlated, at least for those objects 
in which the continuum is unobscured. The PSF magnitude is the closest measurement of the 
nuclear continuum flux. The PSF magnitudes of most low-luminosity AGN are contaminated by 
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their host galaxies. However, if we can find some AGN with minimum host galaxy contamination, 

we can regard their PSF luminosity approximately as their nuclear continuum luminosity, and 
use them to obtain the relationship between emission line luminosity and the nuclear continuum 
luminosity. This relationship can then be applied to other AGN. 

We start from a list of quasars selected by their colors from the SDSS (Richards ct al. 2002) 
(not necessarily in our AGN sample). The fact that they can be selected by their colors means 
that their host galaxy luminosities are weak compared with their nuclear luminosities. We apply 
the stellar subtraction procedure as described in Paper I to these quasars, selecting objects with 
very weak or undetectable stellar absorption line components and use these quasars to establish 
the emission line luminosity vs. nuclear luminosity relationship. 

Figure 2 plots the Ha luminosity vs. the r band PSF absolute magnitude for these quasars. 
The relationship between the two variables is very tight, demonstrating that the Ha luminosities 
for these quasars are indeed excited by their nuclear luminosities. A linear fit to these data gives: 

log(L(Ha)/L0) = -(0.419 ±0.010) - Mr - (0.209 ±0.30) (2) 

The slope is very close to —0.4, which is what we expect in the ideal case, i.e., constant Ha 
equivalent width. The error of the y-intercept does not come directly from the fitting. Instead, we 
directly fit the upper and lower limit of the relationship from the plot (dashed lines in Figure 2) 
and quote 1/2 the difference as the error. If we assume that the SEDs of high luminosity quasars 
and low luminosity Seyfert 1 nuclei are similar, then Equation 2 can be used to infer the r band 
nuclear luminosity for our AGN from their Ha emission line strength. The host galaxy luminosity 
can then be inferred by subtracting the nuclear luminosity from the overall galaxy luminosity. As 
the AGN nuclear luminosity decreases, the overall galaxy fiux decreases accordingly. The minimum 
nuclear luminosity for the AGN to be selected corresponds to the overall galaxy magnitude equal 
to the galaxy target criterion: r-band Petrosian r = 17.77. This places an extra constraint on the 
detection probability obtained from spectral features only. 

Equation 2 is obtained from a sample of quasars, with the measured Ha flux including both 

broad and narrow components. Therefore, the above description of evaluating nuclear luminosity 
only works for broad-line AGN. For narrow-line AGN in our sample, fortunately, the non-stellar 
continuum contribution is negligible compared to the host galaxy luminosity (Kauffmann ct al. 
2003). Therefore, for narrow-line AGN, we simply skip the above analysis and keep the detection 
probability obtained from the spectral features only. 

3. Measuring the Luminosity Function 

The AGN luminosity function $(i>) is defined as the mean number of AGN per unit volume, per 
unit luminosity. In making the plot of the luminosity function without covering several magnitudes 
of luminosity, we also define the auxiliary quantity ^{L): the distribution of AGN per unit volume 
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per log luminosity, given by: 



*(L) = -#4. = ^L*(L), (3) 

d(logio L) logio e 



where ^{L) is the cumulative luminosity function, and is related to by ^{L) = d^{L)/dL. 

$(L) is widely used when making luminosity function plots. Throughout this paper, the luminosity 
functions are all plotted in this form. 

The luminosity function can be measured via various methods. In this paper, we will use 
the non-parametric method which does not require any assumption about the luminosity function 
shape. 

The likelihood that each AGN z, with detection probability function Pi{L) (described in §2) 
have luminosity L, given its redshift and host galaxy properties, is: 

^ pdVmL)dL 
/$(L)p,(L)dL- 

The maximum likelihood solution $(-L) is the function that maximizes the overall likelihood 

$(L) is obtained by minimizing S = — 21n£. In this paper, we use this method to obtain non- 
parametric, linear-interpolated stepwise (Koranyi & Strauss 1997) luminosity functions. The de- 
tailed derivation can be found in Efstathiou, Ellis & Peterson (1988), Koranyi h Strauss (1997) 
and Blanton (2000). 

The luminosity function ^{V) derived above is not normalized. The normalization is done by 
directly measuring the number density of AGN, expressed as: 



n 



where V is the volume spanned by the galaxies in the sample and the sum is over all AGN included 
in this luminosity function calculation. The error associated with the normalization is: 



The final normalized luminosity function is then ^(-L) = n ■ ^{L). 
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4. Luminosity Function Results 

4.1. Overall Ha Luminosity Function 

The overall Seyfert luminosity function can be obtained by combining the Scyfcrt 1 and Seyfert 
2 samples. Here we choose to use the narrow-line AGN sample selected via Kewley's criteria, 
because it is not significantly contaminated by star formation activity. The overall AGN sample 
includes about 4400 AGN: 1317 Seyfert Is and 3074 Seyfert 2s in the redshift range < z < 0.33. 
Since objects at higher redshifts have lower S/N and therefore higher uncertainties in the evaluation 
of the selection function, we limit our measurement of the luminosity function to < ^; < 0.15; this 
excludes about one third of the sample. The selection functions for these Seyferts takes into account 
whether an object is classified either as Seyfert 1 or Seyfert 2 at each luminosity, which includes the 
cases when the classification changes between Seyfert 1 and Seyfert 2 with changed emission line 
strength. Using the methodology discussed above, we obtain the overall Ha luminosity function 
and plot it in Figure 3. 

The Seyfert Ha luminosity function covers a luminosity range of a factor of 10'^. We fit the 
luminosity function with a two-power-law formula as used in the 2dF quasar luminosity function 
by Groom et al. (2004): 

'""WW 

where 0*(L*), L*, a and (3 are free parameters. The best fit results and the of the fit are listed 
in Table 1. We can also fit the curve to the Schechter (1976) function: 

and the best fit (p*{L^), L^,, a are listed in Table 1 as well. Noticing that the luminosity function 
shape is also close to a single power-law, we fit it to a single power-law form: 

Again, the best fit (^*(L*), L*, a are listed in Table 1 (All following luminosity functions are fit 
with the three models and the results are all listed in Table 1). The three best-fit functions are also 
plotted in Figure 3. The shape of the AGN overall Ha luminosity function is closer to a Schechter 
function or a two-power law shape, even though it does not differ much from a single power-law 
(see values in Table 1 for comparison). 

The AGN number density can be estimated by integrating the AGN luminosity function over 
the luminosity range it covers. Using the best-fit Schechter function and within the Ha luminosity 
range of 10^ Lq to 10^ Lq, we obtain the AGN number density as 0.018 Mpc~^. To compare this 
number with the galaxy number density, we integrate the galaxy luminosity function obtained from 
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the SDSS (Blanton et al. 2003b) over an absolute magnitude range of —14 < Mr < —24 and 

obtain the galaxy number density to be 0.094Mpc~^ in the same cosmology. Thus active galaxies 
in this luminosity range comprise about 19% of all galaxies (Ho et al. 1997b; Miller et al. 2003; 
Brinchmann et al. 2004). 

Similarly, we can also obtain the luminosity density by integrating L(j){L) over the entire 
luminosity range. Since the overall Hct luminosity function is a mix of broad Ha plus narrow Ha 
for Seyfert Is and narrow Ha only for Seyfert 2s, and since narrow-line AGN contribute little to 
the r band luminosity, it is appropriate to apply the integration to just the broad-line AGN Ha 
luminosity function. This LF is derived in a similar way as the overall Ha luminosity function, 
and the result is shown in the inserted plot of Figure 3 (thick line with open circles). The Ha 
luminosity function for all Seyferts is over-plotted (thin line with crosses) for comparison. Seyfert 1 
Ha luminosity is very close to the overall Seyfert Ha luminosity even at the low luminosity range. 
The AGN Ha luminosity density integrated over this luminosity function is 1.54 x IO^LqMpc"^. 
Using the relationship between Ha luminosity and AGN r-band continuum luminosity (Equation 2) , 
we obtain the AGN r-band luminosity density as 1.11 x IO^LqMpc"^. Comparing with the galaxy 
luminosity density obtained from Blanton et al. (2003b): 1.84 x IO^LqMpc^^, we find that the 
nuclear activity contributes about 6 x 10~^ of the total light of galaxies in the r-band. 

4.2. Seyfert 1 vs. Seyfert 2 luminosity function 

Our AGN sample is also sufficiently large to evaluate the narrow- and broad-line AGN lumi- 
nosity functions separately. The comparison of the two can be used to demonstrate the relative 
number ratio of Seyfert Is and Seyfert 2s. For this task, the overall Ha luminosity is no longer 
an appropriate parameter for the luminosity function measurements. According to the unifica- 
tion model, the observed overall Ha luminosity in Seyfert Is includes both that emitted from the 
broad-line region and that from the narrow- line region, but the broad component is obscured in 
Seyfert 2s and only the narrow component is observed. Therefore, the same measured overall Ha 
luminosities for broad-line and narrow-line AGN could correspond to different nuclear continuum 
luminosities and therefore the differences of the overall Ha luminosity function of the two types of 
AGN do not necessarily reflect the real differences of the number density of Seyfert Is and Seyfert 
2s (see the inserted plot of Figure 3, Seyfert 1 luminosity function dominates even at the low lu- 
minosity range). A good luminosity function parameter would be emission line luminosities that 
are isotropic to Seyfert Is and Seyfert 2s. In our study, the Ha narrow-line component, [OIII] and 
[Oil] luminosities are good candidates. 

Since we have two narrow-line AGN samples: one selected via Kewley's criteria, which are 
called Seyfert 2 sample; and another one selected via Kauffmann's criterion, which we call Seyfert 
2* sample, we will first compare the Seyfert 1 luminosity function with the Seyfert 2 luminosity 
function, then with the Seyfert 2* luminosity function. In obtaining the luminosity functions, 
the selection functions of objects in each sample are evaluated with the selection criteria only 
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corresponding to this sample in question. Therefore, unhke the case of the total Seyfert sample, the 
detection probability is calculated for Seyfert Is and Seyfert 2s separately. Figure 4, 5 and 6 plot 
the Ha narrow- line component, [OIII] and [Oil] emission- line luminosity functions for Seyfert Is 
(crosses) and Seyfert 2s (open circles) respectively. The differences between the luminosity functions 
of the two types of AGN all show the same pattern: they are comparable at low luminosity, but 
at high luminosities, the Seyfert 2 luminosity function drops off more quickly than does that of 
Seyfert Is. This demonstrates that the number ratio of Seyfert Is and Seyfert 2s is a function of 
luminosity: at low-luminosity, the number density of the two types of AGN are about the same, 
but at high luminosity, Seyfert Is outnumber Seyfert 2s. 

One big concern of our argument of the relative number ratio of Seyfert Is and Seyfert 2s 
comes from the way we define our Seyfert sample. In particular, in Seyfert 2* sample, which is 
selected via the Kauffmann's criterion, we have about 10,700 narrow-line AGN. This is almost 3 
times larger than the Seyfert 2 sample selected via the Kewley's criteria. Our Seyfert 1 sample 
however is not affected by this change of criteria. Therefore selecting AGN via different criteria 
could dramatically change our conclusions about the number ratio of Seyfert Is and Seyfert 2s. 

To test this, we further measure the [OIII], Ha and [Oil] luminosity function for the Seyfert 
2* sample (Figure 4, 5 and 6 open triangles), which is selected via Kauffmann's criterion. Because 
of the large number of objects included in the Seyfert 2* sample, and because our determination of 
the selection function typically requires several hundred Monte-Carlo simulations per object, we are 
unable to measure the probability function for each galaxy as we did for our Seyfert 1 and Seyfert 
2 samples. Instead we randomly select 1/8 of the sample and measure their selection functions, 
being sure to multiply the overall normalization by 8. 

Except for the high luminosity end, the Ha and [Oil] luminosity function of the Seyfert 2* 
sample is larger than that of Seyfert 2s by a factor of about 4. The Seyfert 2* [OIII] luminosity 
function, on the other hand, is only slightly larger than that of Seyfert 2s at intermediate luminosity. 
This is consistent with results of Kauffmann et al. (2003). In figure 2 of the paper, it can be seen 
that most of the Seyfert 2s that are included by Kauffmann's criterion but not by Kewley's criterion 
are those that have intermediate [OIII] luminosities. The results shown in Figure 4, 5 and 6 agree 
with the analysis described in Paper I and Section 1: Seyfert 2*s selected via Kauffmann's criterion 
tend to include all AGN activity, but suffer from contamination from star formation. The only 
exception is the [OIII] emission line, as the [OIII] contamination from star formation is relatively 
weak. Kewley's criterion on the other hand only selects those galaxies that are dominated by AGN 
activity. Therefore, due to the contamination, the Ha and [OH] luminosity functions of the Seyfert 
2* sample include a star formation component and are significantly higher than those of Seyfert 2 
sample, while the [OIII] luminosity function more or less reflects the real AGN distribution. When 
comparing the Seyfert 2* [OIII] luminosity function with that of Seyfert Is, we obtain the same 
conclusion for the relative ratio of Seyfert Is and Seyfert 2s as before: they are comparable at low 
luminosity, but Seyfert Is outnumber Seyfert 2s at high luminosity. 
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A further test for the robustness of our luminosity function results can be applied by constrain- 
ing the Seyfert 1 sample to include only those objects whose narrow Ha component also satisfy the 
Seyfcrt 2 selection criteria as well (see Figure 10 in Paper I). We refer to such objects as Seyfert 
l*s. The probability function for the sample is obtained by requiring the simulated AGN to satisfy 
both the Seyfert 1 and Seyfert 2 criteria. Figure 7 shows the comparison of the [OIII] luminosity 
functions obtained from the Seyfert 1, Seyfert 2 and Seyfert 1* samples. The two broad-line AGN 
sample luminosity functions are almost the same; in particular, the comparison of the Seyfert 2 
and Seyfert 1* luminosity functions shows the same behavior described above. Thus our conclusion 
of the relative number ratio of Seyfert Is and Seyfert 2s is not sensitive to the selection criterion 
uncertainties. 



5. The Relationship between AGN and Host Galaxy Luminosity 

The selection functions we have obtained are based on the assumption that the nuclear lu- 
minosity is independent of the host galaxy luminosity. We can check this assumption by plotting 
the Ha emission line luminosity against the host galaxy Petrosian luminosity for our AGN sample. 
Figure 8(a) shows the plot for Seyfert 2s. There seems to be a strong correlation between the two 
variables: brighter host galaxies appear to have stronger Ha emission lines. But the apparent lack 
of galaxies at the lower-right corner of the plot is due to a selection effect: if the host galaxy is 
too luminous relative to the nuclear emission line strength, the host galaxy will overwhelm the 
active nuclear features and the object cannot be recognized as an AGN. The lack of objects on the 
upper-left corner, on the other hand, simply reflects the fact that AGN nuclear luminosities have 
upper limits: their Eddington luminosities. If we plot these AGN in small redshift ranges over 
which the selection effect is weak, e.g. as shown in Figure 8(b), we see little correlation between 
the two variables. We now carry out a more quantitative test, which involves using our luminosity 
function results. 

Imagine that our assumption is wrong, and for Seyferts, the nuclear luminosity is generally 
stronger in a more luminous host galaxy. Then based on the wrong assumption, the obtained de- 
tection probability function will be underestimated at low luminosity and overestimated at higher 
luminosity. Therefore, the resultant luminosity function will reflect a lower density of low-luminosity 
objects and higher density of high-luminosity objects than there should be. If we use this incorrect 
luminosity function to evaluate the expected luminosity distribution of a subgroup of objects with, 
for example, low redshift or low host galaxy luminosity, the resultant distribution will be system- 
atically shifted to higher luminosity than the observed distribution. Similarly, for a subgroup of 
objects with high redshift or high host galaxy luminosity, the calculated luminosity distribution 
will be shifted to lower luminosity than the observed distribution. Only if our assumption is correct 
will the expected luminosity distribution match up with the observed luminosity distribution for 
any subgroup of our sample. Based on these considerations, we apply the following tests to each 
of our AGN samples. 
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We divide our Seyfert samples (either a Seyfert 1, Seyfert 2 or Seyfert 2* sample) into several 

subgroups by their rcdshifts. For each subgroup, we use our Seyfert luminosity function results to 
evaluate the AGN luminosity distribution wc expect to observe (i.e., the number of Seyferts in the 
samples as a function of redshift), which can be written as: 



where the sum is over all Seyferts in a subsample. We compare FexpectiL) with the observed number 
of Seyferts in this subsample between L and L + AL. Figure 9(a) demonstrates the comparison 
using the Seyfert 2 Ha luminosity function. The observed and expected luminosity distributions 
agree very well with each other. It is not surprising that the two match perfectly for the entire 
Seyfert 2 sample, since the luminosity function is measured from this sample. But the fact that 
they agree in each redshift bin demonstrates that our luminosity function reflects the true Seyfert 
2 distribution. In addition, this demonstrates that our assumption of the independence of AGN 
and host galaxy luminosity is reasonable. 

We next run a similar test by dividing the Seyfert 2 sample via host galaxy luminosity. By 
doing so, we can directly examine the relationship between nuclear luminosities and host galaxy 
luminosities. Figure 9(b) shows the result of this test for the Seyfert 2 Hq luminosity function. 
Higher luminosity Seyfert 2 galaxies have more luminous Ha emission lines (as they tend to be at 
higher redshift on average), as we saw in Figure 8(a). However, in each subgroup, the observed 
and expected luminosity distributions match up impressively well, which gives us a quantitative 
confirmation that the nuclear luminosities for Seyfert 2s are indeed not strongly correlated with 
the host galaxy luminosities. 

In Figure 10, we show the tests done for the Seyfert 1 sample via the [OIII] luminosity function. 
It is still a perfect match between the observed and expected luminosity distribution for every red- 
shift and host galaxy magnitude subgroup. This demonstrates that the host and AGN luminosities 
are not strongly correlated for Seyfert Is either. 

We apply similar tests for other luminosity functions, and always found good match between 
the expected and observed luminosity distribution. 

As a further test of the correlation between host galaxy and nuclear luminosity, we apply 
a correlation test to our Seyfert samples. Correlation tests for truncated data sets have been 
developed by Efron Sz Petrosian (1992) and Maloney & Petrosian (1999). Our sample is a one-side 
truncated sample in flux. Fan et al. (2001) discussed the correlation test for a one-side truncated 
sample with non-sharp selection probability, which is appropriate in our study. 

For each of our AGN samples, we have a one-side-truncated data set {Ln^i, L^^i}, where Ln,i 
and L/j j are the emission line luminosity and host galaxy luminosity of the ith AGN in the sample. 
For each AGN in the sample, we can define a comparable data set as: 




(11) 



Ji — {j '■ Lin,j > Ln,i} ■ 



(12) 
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We denote the number of points in the set as Ni. Each AGN in the data set has a detection 
probabihty p{Ln^i, L^^i), thus the total number of AGN that have nuclear luminosity larger than 
Ln,i, taking into account those that are missed by selection effects, should be: 

rp. Pj{^n,i-iLhj) ^^^^ 

Assuming that L„ and are independent, we sort Tj by their host galaxy luminosities. The rank 
Ri of L/i j among the Tj objects, which is defined as: 

should be distributed uniformly between and Ti. The expectation value of the distribution Ri is 
then Ei = Ti/2 and its variance Vi = Tf /12. We can construct Kendall's r-statistic: 

r = ^^^IZ^^^ (15) 

If |r| < 1, the nuclear luminosity and host galaxy luminosity are not correlated at the la level and 
can be treated as independent. 

We assume a relationship between the host galaxy luminosity and nuclear luminosity with a 
power-law index a: 

Lh = K (16) 

and run the {Lh,Lf^} data set through the correlation test to find the value for a (when r = 0) 
and its error (the differences of a when r = and |r| = 1 ). Table 2 lists the results for Seyfert Is, 
Seyfert 2s and the full Seyfert sample when Lnudear is represented by [QUI] and Hcc luminosities 
respectively, a 1 in all cases. We conclude that the nuclear luminosity is essentially independent 
of the host galaxy luminosity for all Seyferts in our sample. 

The above statement is obtained from our AGN sample, which unfortunately covers only a 
small range of host galaxy luminosity. Outside of this luminosity range, the independence of the 
host and nuclear luminosity might not hold. For example, Kauffmann et al. (2003) suggested that 
AGN of all luminosities reside almost exclusively in massive galaxies. Within this regime, there is 
only a very weak dependency of nuclear luminosity and galaxy mass, which is in agreement with 
our result. But the probability of a galaxy to host an AGN drops dramatically at galaxy masses 
lower than a few xlO^^ Mq. 



6. Comparing LFs with Previous Work 



Most previous AGN luminosity function results obtained in the literature have been functions 
of B band magnitude. To make comparisons of our luminosity function results with previous 
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investigators, we need to relate the emission line luminosity with the AGN nuclear continuum in 
the B band. In §2, we have found that quasars with little host galaxy contamination show a tight 
relation between their emission line luminosities and r band PSF magnitudes. We assume that all 
AGN follow the equivalent relationship found with the i-band absolute magnitude: 

log(L(Ha)/Lo) = -(0.470 ± 0.011) • Mj - (1.38 ± 0.41) (17) 

Again, as in equation (2), the error in the y- intercept is chosen to include the la scatter of the 
relationship. With this equation, Ha luminosity function can be converted to the i band luminosity 
function. Based on the general quasar color model (Schneider et al. 2002, 2003; Hopkins et al. 
2004), we found g — i ^ 0.55 at around z ~ 0.1. Also noticing that for normal quasars, B g, we 
can use the simple relationship: 

5 ~ i + 0.55 (18) 

to convert the i band luminosity function further to B band. 

Figure 11 shows the converted B band magnitude luminosity function for all Seyferts. The 
uncertainties in equation (17) arc also shown (dashed line). The AGN luminosity function results 
of Huchra & Burg (1992) are plotted as triangles, while those of Ulvestad & Ho (2001) are plotted 
as stars. Both results have been converted to Ho = 100km Mpc~^. Their choice of {Q,m,^A) 
are also different from ours, but this won't make much difference at our redshift. In the region 
of overlap, these results are in good agreement with ours, although ours reach to much fainter 
luminosities. Wc should bear in mind that Huchra & Burg and Ulvestad & Ho used host galaxy 
luminosities without isolating the nuclear luminosity as we did. Therefore their luminosity function 
is much flatter than ours since low-luminosity AGN can easily be buried in host galaxy luminosities 
and be undetected. But since their sample basically covers the high-luminosity end of the AGN 
population, their luminosity function is not severely contaminated. 

To compare our luminosity function results with high redshift quasar luminosity functions, 
we take the 2dF QSO luminosity function obtained by Groom et al. (2004) and converting it to 
our assumed cosmology. In Figure 12, we plot our B band AGN luminosity function (solid line 
with dotted lines as error range) and converted 2dF QSO luminosity function (dashed line with 
triangles). Our low-rcdshift AGN luminosity function agrees very well with the quasar luminosity 
function for the lowest redshift range. However, our data show that the faint end of the 2dF QSO 
luminosity function at low redshifts does not flatten out. This reveals that the color selection of 
quasars is subject to incompleteness at the low-luminosity range where the host galaxy dominates 
(Boyle et al. 2000). Unfortunately, the large uncertainty on our luminosity function, due to the 
scatter of the conversion from the Ha luminosity to the B band magnitude, makes it difficult to 
draw conclusions about AGN evolution from 2; ~ 0.1 (our result) to z 0.4 (the lowest redshift 
of 2dF QSO LF). The ongoing effort of getting the quasar luminosity function from the SDSS will 
test this directly. 



7. Reddening Correction 



The luminosities measured in this paper are not corrected for any reddening of the AGN 
by material in the vicinity of the AGN or the host galaxy. However, reddening is an important 
factor that can dramatically change the observed luminosity of an AGN. The standard method of 
reddening correction is to measure the Balmcr decrements and assume an intrinsic Ha/H/? value 
(e.g. 2.86) and a reddening function (e.g. Scaton 1979; Chariot & Fall 2000). However, AGN with 
large Balmer decrements have been observed before (Anderson 1970; Baldwin 1975; Weedman 1977; 
Rix et al. 1990; Goodrich 1990; Barcons et al. 2003), many of these can not be simply explained 
by dust obscuration. Instead, AGN might have a higher intrinsic Ha/H/? ratio, probably because 
of higher densities and radiative transfer effects. In addition, the relative distribution of dust and 
emission line gas can be quite complicated, while most models assume that the dust lies in a screen 
in front of the line-emitting regions. Due to these difficulties and uncertainties, we believe it is not 
appropriate to do a standard reddening correction using observed Balmer decrements, even though 
reddening is clearly present. The whole issue of AGN obscuration needs much further investigation 
and multiwavelength observations from infrared and X-ray, and is beyond the scope of this paper. 
We will leave the luminosity function derived in this paper reddening uncorrected until reliable and 
accurate AGN reddening correction methods are developed. 



8. Summary and Discussion 

Using the AGN sample obtained from Paper I, we have evaluated the AGN [OIII] and Ha 
emission line luminosity functions. Because of the complexity of our AGN selection criteria, we use 
a Monte-Carlo simulation for each AGN to obtain its selection probability function. The overall Ha 
emission line luminosity function for Seyferts covers a broad range of luminosity and approximately 
follows a two power-law form. After converting to B band magnitude luminosity functions, we found 
that our results are comparable to the luminosity function obtained by Huchra et al. (1992) from 
the CfA redshift survey and Ulvestad & Ho (2001) from the Palomar spectroscopic survey of nearby 
galaxies. It also agrees well with the low-redshift bin of the 2dF QSO luminosity function by Groom 
et al. (2004). By integrating the overall Ha luminosity function from 10^ to 10^ Lq, we obtain the 
AGN number density of 0.018 Mpc~^, approximately 20% of all galaxies. We also estimated the 
AGN r band luminosity density as 1.11 x lO^LoMpc"^, which is about 6 x 10~^ of the luminosity 
density of galaxies. 

The comparison of the [GUI] luminosity functions for Seyfert Is and Seyfert 2s reveals that 
the relative number ratio of the two types of AGN is a function of luminosity: at low luminosity, 
the ratio of Seyfert Is to Seyfert 2s is about one, while at high luminosities, Seyfert Is outnumber 
Seyfert 2s. The conclusion is not sensitive to the uncertainties in the selection criteria. The 
paucity of narrow-line objects in high luminosity AGN is unexpected in the simplest unification 
model, and various models have been proposed to accommodate it. For example, Lawrence (1991) 
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proposed that the opening angle of the dust torus is larger for more luminous AGN, presumably 
because increased luminosity means more dust sublimation, and thus the broad-line region can be 
seen over a larger opening angle. Konigl & Kartjel (1994) proposed a disk-wind AGN model, in 
which the radiation pressure force flattens the dust distribution in objects with comparatively high 
bolometric luminosities, thus the opening angle of the dust torus becomes larger for more luminous 
AGN. Our result for the first time gives the opportunity to determine the detailed opening angle 
as a function of luminosity and will provide detailed constraints to various models. This will be 
fully investigated in future papers. To fully understand the physics behind the lack of Seyfert 2s 
at high luminosity, multiwavelcngth observations of AGN, especially the infrared will be necessary. 
By directly measuring dust emission in Seyfert Is and Seyfert 2s, various dust properties can be 
determined, and will help to understand if dust is playing a major role in the difference between 
Seyfert Is and Seyfert 2s. 

Our AGN luminosity function is obtained by assuming that the AGN nuclear luminosities are 
independent of the host galaxy luminosities. Detailed verification of this assumption has been car- 
ried out and we show strong evidence that the two variables are not strongly correlated. AssTiniing 
that AGN are hosted by normal galaxies, their host galaxy bulge luminosities are well correlated 
with their velocity dispersions, which in turn are strongly correlated with the central black hole 
masses (Gebhardt et al. 2000; Ferrarese &; Merritt 2000; Ferrarese et al. 2001; Tremaine et al. 
2002) . The independence of the nuclear luminosity and the host galaxy luminosity implies that the 
nuclear luminosity does not depend on the central black hole mass (Heckman et al. 2004). There- 
fore, for the AGN in our sample, the principal variable that determines the nuclear luminosity is 
the Eddington ratio, defined as A = Lhoi/LEddington, instead of the black hole mass. 

The independence of AGN nuclear luminosities and their central black hole masses is sup- 
ported by several observations. For example. Woo & Urry (2002) compiled a sample of 377 AGN 
with known black hole masses and calculated bolometric luminosities for 234 AGN among them. 
When relating these bolometric luminosities with their black hole masses, they found no significant 
correlation between the two variables. O'Dowd &; Urry (2002) also reported a very weak trend 
between the nuclear and host luminosities for a sample of radio loud AGN. 

However, the same statement does not seem to work for high-redshift quasars. A black hole 
mass inferred from the observed luminosity is a lower limit. Thus given their high luminosities, 
one would infer implausibly high black hole mass (> IO^^Mq) if the most luminous quasars are 

accreting at substantially sub-Eddington rates, (see Vestergaard et al. 2002 for black hole mass 
measurements for high-redshift quasars). Thus, it will be interesting to explore the transition from 
low to high Eddington rate accretion (McLure &: Dunlop 2004). Hopefully, the ongoing evaluation 
of the quasar luminosity function obtained from the SDSS will yield insights on this issue. 

Our result will also bring insights to X-ray background and UV background problems. Detailed 
investigation is beyond the scope of the current paper. The multiwavelength observation of the 
sample will be explored in future papers. 
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Fig. 1. — Three examples of our selection function measurements. The left column shows the 
original spectra of three AGNs and their simulated spectra after the emission lines have been 
scaled down by a factor of 0.4 (offset by -10). The right column plots the corresponding detection 
probability functions measured. 
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Fig. 2. — The r band PSF magnitudes vs. the Ha luminosities for quasars with Httle or no 
stellar contamination. For these quasars, the PSF magnitudes can be approximately regarded as 

all coming from the nuclei. Therefore, the relationship reflected in this plot can be regarded as 
the relationship between the nuclear luminosity and the Ha emission line luminosity. The solid 
line is Equation 2. The dashed lines indicate the approximate scatter of the relationship and the 
differences between dashed lines and the solid line are considered as the error of Equation 2. 
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Fig. 3. — The overall Ha luminosity function for Seyferts, including both Seyfert Is and Seyfert 
2s. We fit the luminosity function with Schechter function (Equation 9), Two- Power-Law function 
(Equation 8) and One-Power-Law function (Equation 10). The best fitting parameter and corre- 
sponding are listed in Table 1. The inserted plot shows the Seyfert 1 Ha luminosity function 
(thick line with open circles) compared with the overall Ha luminosity function. 
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Fig. 4. — The narrow-line component Ha luminosity functions for Seyfert 2s (open circles, selected 
via Kewley's criterion), Seyfert Is (crosses) and Seyfert 2*s (triangles, selected via Kauffmann's 
criterion) separately. The Seyfert 1 and Seyfert 2 luminosity functions are about the same at low 
luminosity. But at high luminosity, the Seyfert 2 luminosity function drops off more quickly than 
does the Seyfert 1 luminosity function. The Seyfert 2*s luminosity function is significantly larger 
than that of Seyfert 2s. To avoid crowding the plot, we didn't show errorbars of the Seyfert 2* 
luminosity function, which are comparable to the errorbars of the other two luminosity functions. 
The same applies to the following three plots. 



-26- 




105 IQ6 IQ7 IQB 

[OIII] Luminosity (Lq) 



Fig. 5. — The [OIII] luminosity functions for Seyfert 2s (open circles), Seyfert Is (crosses) and 
Seyfert 2*s (triangles) separately. The comparison of the Seyfert 1 and Seyfert 2 luminosity func- 
tions shows us that the number ratio of Seyfert Is and Seyfert 2s is a function of luminosity. At low 

luminosity, the number density of Seyfert Is and Seyfert 2s arc about the same, but at high lumi- 
nosity, Seyfert Is gradually outnumber Seyfert 2s. The fact that the Seyfert 2* luminosity function 
is comparable to that of Seyfert 2s demonstrates that we are not missing much AGN activity in 
the Seyfert 2s sample. The comparison of the three confirms our conclusion of the relative number 
ratio of Seyfert Is and Seyfert 2s. 
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Fig. 6. — The [Oil] luminosity functions for Seyfert 2s (open circles), Seyfert Is (crosses) and 
Seyfcrt 2*s (triangles) separately. The comparison of the Seyfert 1 and 2 luminosity functions 
shows the same pattern as of Figure 5. Again, the Seyfert 2* luminosity function is significantly 
larger than that of Seyfert 2s because of star formation contamination. 
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Fig. 7. — The [OIII] luminosity functions for Scyfcrt 2s (open circle), Scyfcrt Is (crosses) and 
Seyfert Is (triangles) whose narrow components also satisfy the narrow-hne AGN selection criteria 
(noted as Seyfert Is*). There is little differences between the luminosity functions for Seyfert Is 
and Seyfert Is* and the comparison of the luminosity functions of Seyfert 2s and Seyfert Is* reflects 
the number ratio of the two types of AGNs free of the selection criteria uncertainties. 



-29- 




>#< 



-16 -18 -20 -22 -24 -16 -18 -20 -22 -; 

"^petrosian "^petrosian 

Fig. 8. — The r band Petrosian magnitudes vs. the Ha luminosities for the narrow-hne AGN 
sample. It appears that there is a correlation between the two luminosities. But the apparent 

correlation is very likely to be due to selection effects. When we choose a subsample with small 
rcdshift range, so that the selection effects are not significant (shown in the right figure), we find 
little correlation between the two luminosities. 
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Fig. 9. — The observed luminosity distribution for Seyfert 2s (liistogram) compared with the 
expected luminosity distribution (solid line) calculated from the Seyfert 2 Ha luminosity function 
(Figure 4) for different subgroups of the Seyfert 2s sample. Left: the entire Seyfert 2 sample 
(shown in the top panel) is grouped by their rcdshifts. Right: the entire Seyfert 2s sample is 
grouped by their absolute r band Petrosian magnitudes. The match up between the observed 
luminosity distribution and the expected luminosity distribution in each subgroup demonstrates 
that our luminosity function results are correct and the assumption that the host galaxy luminosity 
is independent of the nuclear luminosity is reasonable. 
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Fig. 10. — The same test as in Figure 9 for the Seyfert 1 [OIII] luminosity function. The match 
up between the observed luminosity distribution and the expected luminosity distribution in every 
subgroup demonstrates that the independence of nuclear luminosity and host galaxy luminosity 
works for broad-line AGNs as well. 
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Fig. 11. — Our converted B band magnitude luminosity functions for Scyfcrts (black) with errors 
(dashed lines) compared to the AGN luminosity function obtained by Huchra et al. (1992) from 
the CfA redshift survey (solid triangles) and by Ulvestad &: Ho (2001) from the Palomar Seyfert 
Galaxies (stars). Our luminosity function basically agrees with the other two luminosity functions. 
But we should bear in mind that Huchra et al. (1992) and Ulvestad &; Ho (2001) used B band 
magnitudes for the whole galaxies in their luminosity function calculation, rather than the nuclear 
B band magnitudes we have used. 
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Fig. 12. — Our converted B band magnitude luminosity function, compared with the 2dF quasar 
luminosity function, converted to Hq = 100km s~^ Mpc~^. From left to right, the different 2dF 
quasar luminosity functions are for the redshift range: 0.40 < z < 0.68, 0.68 < z < 0.97, 0.97 < 
z < 1.25, 1.25 < z < 1.53, 1.53 < z < 1.81 and 1.81 < z < 2.10 respectively. At luminosities fainter 
than Mb = —21.5, the 2dF quasar luminosity function might suffer significant incompleteness due 
to the effect of the host galaxy, therefore the comparison should be made only for luminosities 
brighter than Mb = -21.5. 
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Table 2. Correlation Test Results 
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